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Abstract 



A large set of complex path solutions for the Hartle Hawking semi-classical 
wave function are found for an inflationary universe in the "slow roll" regime. 
The implication of these for the semi-classical evolution of the universe is also 
studied. 

O ! Hartle and Hawking Q suggested that one could overcome the problem of supplying 

initial conditions for the universe by supplying a theory of those initial conditions within 
the context of quantum gravity in the path integral formalism. One of the difficulties with 
any quantum or classical system is that the theory only predicts the evolution of a system- 
ie the correlations between the states of the physical attributes of a system at different times. 
Unfortunately in studying the universe as a whole, there exists no way of determining what 



£jrj| the early state of the universe other than by observations made now. There is no way of 
testing whether or not one's theory of the early universe is actually correct, only of displaying 
what initial conditions must have been present at that early time to produce the universe 
as it is now. Since any physical theory always has a set of initial conditions which can 
produce any arbitrary final state, one is left with aesthetic judgments as to whether or not 
the required initial conditions are pleasing or not. 

However, within the context of the path integral formulation of quantum Einstein gravity, 
Hartle and Hawking demonstrated that one could generate a natural theory of a replacement 
for those initial conditions. Namely, since Einstein's theory allows changes in the structure 
of space-time, one could imagine a requirement on the universe that it not have any initial 
conditions. The universe can simply run out of time into the past, with some unique structure 
demanded at the time that the universe came into existence. 

This suggestion is not tenable within the context of ordinary General Relativity because 
of the various singularity theorems for Lorentzian space-times. However, for a Euclidean 
space-time, Einstein's equations allow singularity free solutions with cross sections repre- 
senting various possibilities for the universe at some instant in time. Such solutions, while 
of dubious relevance for the classical theory, can come into their own in the quantum theory, 
especially the quantum theory formulated in terms of path integrals. 

The Hartle Hawking prescription was thus that one could calculate a complex valued 
function, called the wave-function of the universe, within the path integral formalism. This 
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function is a function of three geometries- interpreted as the geometry of the universe at 
some particular time- together with configurations of classical fields on that that three 
geometry. The definition they suggested was that this wave function was to be calculated 
by means of a path integral in which one integrated over all non-singular four geometries 
which had only this three geometry as a boundary, together with fields on that 4 geometry 
such that the value of those fields on the specified three geometry were the given fields. I.e., 

y(g^), <5>(xeg^) = J e iS ^^ G ^54>8g (1) 

where the four-geometry and the fields <fi defined on that four-geometry, are non- 

singular. 

While this formulation of the wave-function is certainly appealing intuitively, it suffers 
from numerous technical problems. Since path integrals even in standard situations are 
typically taken over highly pathological paths (eg for a free particle over paths which are 
continuous but not differentiable), the requirement that Q^> and (f) be non-singular is difficult 
to specify exactly. Furthermore, this path integral cannot be evaluated even in very simple 
cases. Finally, once one has calculated this wave-function, it is unclear how it should be 
used to actually derive predictions from the theory. 

However the purpose of this paper is not to address these difficulties. Rather we will 
take a standard naive approach to the problem. We will assume that the wave-function 
is well approximated by the semi-classical approximation, in which one assumes that the 
extrema of the action, S will dominate the path integral. Furthermore, we will neglect the 
"determinant" terms, which arise from taking into account the integral over paths which 
are "near enough" to the classical path that one can make the quadratic approximation to 
the action of these paths near the semi-classical path. Finally, we will interpret this wave 
function as giving us a semi-classical Hamilton Jacobi function, which we can then use to 
predict the semi-classical evolution of the universe. 

One's first reaction is that surely this semi-classical evolution will be simply some clas- 
sical solution to Einstein's equation. The difficulty arises in that there exist no classical 
Lorentzian solutions to Einstein's equations which are non-singular and have only the given 
three-geometry plus matter fields as a boundary. The proposal made by HH was there- 
fore that these semi-classical solutions be taken to be a combination of Lorentzian and 
Euclidean four geometries, matched together at some time. The fiducial example is if we 
take the simplest situation, namely one in which the only source for the gravitational field 
is a cosmological constant, and that the geometry be restricted to a maximally symmetric 
geometry, with three geometries being metric three-spheres. In this case one can join a 
De-Sitter Lorentzian universe to a four dimensional metric sphere (a Euclidean solution to 
the Einstein equations with cosmological constant). In order that Einstein equations also be 
valid across the junction between the solutions it is necessary that the extrinsic curvatures of 
the three geometries at the join be continuous. However, given some definition of the time, 
one can characterize the Euclidean solution as having imaginary time, while the Lorentzian 
has real time. The extrinsic curvature in the Euclidean space is therefore imaginary while it 
is real in the Lorentzian, and thus matching can be smooth only if the extrinsic curvature is 
zero [Q]. Thus the join between the two solutions must take place at the minimum 3- volume 
for the De-Sitter space, and the maximum 3-volume for the Euclidean- ie one must match 
the center of the De-Sitter throat to the equator of the Euclidean sphere. 
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This four geometry has the feature that at any "time" the three geometry obtained 
by taking an equal time cross section is a real three geometry, and the vague idea has 
developed that this is typically what one wants to have happen. However, one begins 
to quickly realize that, in general (ie in more complicated situations), this is impossible. 
Essentially, the requirement that one have a real three-geometry at each and every time 
would require a matching between the Euclidean and Lorentzian regimes such that the 
momenta corresponding to each of the fields, and to the geometry would all be zero at the 
junction. If not, the momenta would have to be complex in one or the other of the regimes, 
and would thus quickly (because of the equations of motion) create complex field strengths 
and complex metrics. However, as we shall see in very simple cases, there exist no solutions 
to the classical Euclidean Einstein equations which are both regular and have one three 
boundary on which all of the momenta of all of the fields and of the geometry are zero. 

One is therefore driven to investigating complex three-geometries, and complex four 
geometries as classical solutions to Einstein's equations. (We also advise the reader to 
read the papers by Halliwell and Louko |4j] and Hartle and Halliwell || which contain a 
deeper discussion of the motivation behind finding complex solutions. Note also the paper 
by G. Lyons || who qualitatively examined the case where the potential for is quadratic 
and arrives at some of the same conclusions we do). In this paper we will investigate the 
complex solutions for a simple problem, namely one in which most of the degrees of freedom 
of the geometry and the matter fields are frozen, so that the geometry has a high degree of 
symmetry (in particular having the symmetry of the metric three-sphere) but in which there 
exists a matter field given by a scalar field with a non-trivial potential, but with the same 
spherical symmetry of the metric. The potential for the scalar field will be assumed to be 
almost flat (ie with small slope) and we will examine the complex solutions as a power series 
in that slope of the potential. With zero slope this problem of course simply reduces to that 
originally studied by Hartle and Hawking, since a flat scalar field potential is equivalent to 
a cosmological constant. 

We will find that there are in fact an infinite number of complex semi-classical solutions. 
Many of them will be analytically related to each other. They will also give an action 
for the final state of the system which is the same for large classes of the paths. There 
are however an infinite number of analytically inequivalent paths, many of which lead to 
different wave-functions. 

Furthermore, if we examine the classical solutions implied by these solutions, we find 
that the classical space-time need not be the same as (in fact because of the complex nature 
of the extremal solutions, cannot be the same as) any of the extremal paths. We suggest 
that one interpretation of the Hartle Hawking wave function is that the classical universe 
has no initial singularity because it bounces. Ie, we suggest that the Hartle Hawking wave 
function is a way of imposing the requirement on the quantum system that the universe 
bounce ( perhaps due to quantum corrections) before it hits a singularity. 

I. THE MODEL 

The universe is assumed to have the geometry 

ds 2 = N(t) 2 dt 2 - a(tf (dr 2 + sin 2 (r)(d6 2 + sin 2 (9) dtp 2 )) (2) 
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with a homogeneous scalar field (f>(t). Since we will be interested in complex extrema to the 
action, we will assume that all of N the lapse function, a the scale factor, and the scalar 
field can be complex. However, we will assume that the real universe has real scale factor 
a and real scalar field 0. Thus the wave function of the universe is a function only of real 
a — cif and real = 0/. Ie, we will calculate \&(a/, 0/) via the semi-classical approximation. 
The action will be assumed to be the standard Einstein Hamiltonian action, 



where 



S = J (V J ' 7ij + tt - NH - NiH 1 ) d 3 xdt (3) 
H = (jf*^ (7*7*1 - -V^R+ {\** 2 + V7( 2^=M"7 y + W)) ( 4 ) 

Hi = *J + 1%** + - o,v j (5) 

Because of the symmetry we have assumed for our fields, ifj is automatically equal to zero, 
and we can write the action in terms or a, and their conjugate momenta ix a and ix^ so that 

S = I ir a a + vr^0 - NHdt (6) 



H = - 6a + + a 3 V(0)) (7) 

The Hartle Hawking procedure is often characterized as allowing time to become imag- 
inary. However, we will assume that t is real throughout. The transition to Euclidean 
space-time will be done by allowing N the lapse function to become purely imaginary. How- 
ever, in addition to allowing a and to become complex (as they must) we will also allow N 
to be not only purely real or purely imaginary but also to be an arbitrary complex function 
of the real time t. Ie, we will be looking at complex four-metrics and complex fields on a 
real four dimensional manifold characterized by coordinates t r 9 (p. 

Despite the simplicity of this action, the extrema are still difficult to find analytically. 
We will therefore find them by assuming that V((f>) is almost constant, and can be written 
as 

V{4>) = V + e0 (8) 

We will then assume that the solution is analytic in e and solve our problem as a power 
series in e. Ie, we assume that V(0) is essentially that for inflation, with the solutions we 
examine being for conditions such that a significant era of inflation takes place. 

It will be useful for future work to note that this action has a symmetry. In particular 
if we allow V to go to V + eQ and to go to — Q the action is left invariant. 

The extrema are defined by the equations 

1 da 



N dt 12a 
1 d0_ 1 

1 dix n 



N dt 



(9) 
(10) 

(11) 



< ~ 6 - A*J + 3aV (0) 



24a 2 a 2a 4 



1^1- _ a 3^ (12 ) 

N dt~ a d4> [iZ) 
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We also have the fifth equation, an integral of the motion of the above equations and the 
constraint equation corresponding to the variation of the action with respect to N, of 



H 



24a a v 2a 3 
Using the equations of motion, this can be written as 



(13) 



-6aa' 2 -6a+-a 3 (p' 2 + a 3 V((p) =0 (14) 
2 

where ' denotes j^^. We note that N occurs in these equations only in the combination 
Ndt. We can thus define a new variable r by 



J N(t)dt 



(15) 



where we will take the zero of r to occur at the "time" when a is zero. Since we are allowing 
N to be complex, r will also be a complex function of t. r(t) will be a parameterized 
path through the complex r plane. Note that these different paths are not simply different 
coordinates. A coordinate transformation on t cannot create a change in for example the 
ratio of the real to imaginary part of r. The paths through r will be coordinate invariant 
paths. ' will then denote the derivative with respect to r. 

We want the four geometry to be non-singular at r = 0. Since r = is defined to be 
the initial point along the path where a = 0, the regularity condition on a will be one on 
the derivative of a and on the momenta 7f a and n^. Examining the metric, we see that it 
is regular at a = if at that point a goes as r and r is purely imaginary so as to create a 
Euclidean metric. (For example, the curvature, Gq = 3(a' 2 + I) /a 2 , is non-singular as a goes 
to zero only if a' 2 + 1 goes to zero at least as fast as a 2 does) Thus, the boundary condition 
at r = is that a' = ^ = ±i at r = 0. Furthermore, the boundary condition of ir^ is 
that 7T0 must go to zero, or the equation for 7T a cannot drive it to zero sufficiently quickly at 
r = 0. We thus have the conditions 

a(0)=0 (16) 
a'(0)= pi => ix a = -12t as r -> (17) 

M°)= ( 18 ) 

where f3 — ±1. We also have the final conditions 

«( r /) = a .f = ±M ( 19 ) 
4>{r f ) = 4>f (20) 

where however r/ is a free variable which can be chosen so as to allow these equations to be 
obeyed. As expected we thus have five boundary conditions on the two variables obeying 
second order equations, a, <f> and the third variable r/. Both signs for the final radius a/ are 
valid since it is only a 2 which enters the metric. 
The action in terms of the r variable becomes 

S= f Tf n a a' + 7T 0' - Hdr. (21) 

J 
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We will solve the classical equations of motion in the r variable and as a power series in the 
slope of the linear potential, e = dV j^ ■ At each order in e, the point Tf at which a(rj) = a/ 
will also change. We will define To as the value of Tf for e = 0. If is one of our dynamic 
variables, a or and if qf^ is the zeroth order solution, and 



% = qf ] + E -fq* (22) 
n ! 



Then 



gives 



qi (r f (e)) = q lf (23) 



^(roH-^gfVo) (24) 
S\(t )= -S 2 r f qf\r ) - {5r f fqf\r ) - 25r f 5q'(r ) (25) 

We shall use these equations for qi = a to calculate the expansion of 77. 
Let us now examine the expansion of the action, 

S= f T, (7r t q' t -H)dr (26) 
Jo 

To zeroth order we have 

S<°) = I Trf a^'dr (27) 

We also have, for arbitrary e, 

~d7 = ~d7 fafa^' e ) % ( r /( e )' e )' ~ ^fa' e )] 

+ P + _ (^^£ + ^% + d _E\ dT (28) 

Jo de 1 1 de \ diii de dqi de de I 



f[,w e ), e »„ We )i-(,|);'-('f* <»> 



where we have used the Hamiltonian equations of motion and the constraint equation H = 
to simplify the equations. But, 

&(rM,H*M^ = % = (30) 

so we finally have 

dS(a f ,(f)f,e) rs dHfa, q h e) 



de 

In our case, we have that 

dH 



m cm {7ij,qi,e) 
~ Jo de [6L) 



de (32) 



But, we also have from the equation of motion for <p that 



a 3 = _i (VyY (33) 



Thus 



dS__ 1 

de e 

1 1 r T f 



a 3 0')'0^ (34) 



aj(j)(T f y(f)f - - / a 3 (0') 2 dr (35) 

-0/ f Tf {a 3 <pydr - - r a\(j)') 2 dT (36) 
e Jo e Jo 

, Tf a 3 dr - - r a 3 (0') 2 dr (37) 
e jo 

— dr-- / a 3 (0') 2 ^ (38) 
/o oV Q e Jo 

However, by exactly the same line or argument which led to eqn [H], we have 

dS(a f , (j) f , V , e) rs dH 



-dr (39) 



Thus we finally have 



dV Jo dV Q 



dS dS 1 rs 3 ,jj, 2 j ,An\ 



We can solve this in part by requiring that Vq occur in S only in the combination Vo + e0/. 

For our particular problem, the second term on the right is of order e. Thus, the first 
order term SS is determined by the zero order term, and the first interesting terms is the 
second order term 5 2 S. Since <$' = e5<f)' + 0(e 2 ), the second order term in S is determined 
by the zeroth order term in a and the first order term in 5(f). 

We will look at the first order solutions for both a and 0, namely 8a and 5(f). 



II. SOLUTIONS AND ACTIONS 
A. Zeroth Order solutions 

With our potential, V(<p, e) = Vo + £(f> we have the equations of motion 

- 6aa' 2 - 6a + ^a 3 /2 + a 3 (V + e<j))= (41) 

(a 3 (f/y= -a 3 e (42) 

The zeroth order solutions give that a 3 (f>' a constant. If this constant were non-zero, then 
(f> would diverge at r = at least as 1/r 2 and the solution would not be regular. Thus, to 
zeroth order, is a constant, and from the final boundary condition <p(r) = (f>f. 
The zeroth order solution for a of 
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a (0)« = a (o)2 (H/6) _ ! (43) 
is (with the appropriate boundary condition) 

a (0) (r)=/9^sin(^r) (44) 

We will define v = since it will occur time and again, so the solution becomes 

a (°)( r ) = 0- sin(w/T) (45) 

This solution obeys the boundary condition on a at r = 0. We also have the boundary 
condition a^(ro) = a/. We readily see that this equation has a number of solutions, since 
sm(ivr) is periodic in imaginary r with period 2ni/i>. Since a/ is real, we have, defining 

r R = Real(r Q ) (46) 
77= Imag(r ) 



that r/j and 77 obey 



/3 cosh (vtr) sin (^r/) =va$ (47) 
cos(i/t/) = (48) 



For cif > 1/u, this has solutions for tj 

Tr = (2n + sign(a f )P-)- (49) 

for all n. In addition, the equation for t r 

cosh(z/r fi ) = v\df\ (50) 

has both positive or negative solutions for tr. In what follows an important variable will be 
cos(iz/T ) whose sign depends on both the sign of ay and on the sign of t r . I will define the 
sign, 7, (ie, 7 = ±1) by 

cos(iur ) = 7 [+i\^Jv 2 a 2 f - 1|) (51) 

for z/a/ > 1. In this case 

7 = sign(r i? )sign(a / ) (52) 
For vdf < 1, tr must be zero, and tj, the imaginary part of r must obey 

sin(z/r/) = —j3uaf (53) 

and 7 is given by 
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cos(iut ) = cos(vtj) = 7 (+|y 1 — z/2a /l) • (54) 

In most of what follows we will assume that we are dealing with uaf > 1. 

Note that for tr = and 77 = mir/v, a goes to zero, m = corresponds to the initial 
value r = 0. For even m, these zeros will be repetitions of the initial state. However, for 
odd m these zeros will turn out to be singularities in the solutions of the equations at higher 
orders of e. 
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Figure 1 

The location of the zeros of a(r) and various possible endpoint values for a given a/ of the 
paths in complex r space. 

Figure 1 shows these poles, the points where a is zero and values of r which satisfy the 
l a ^( T o)l = \ a .f\- i n the complex r plane. This figure also also indicates a few possible paths 
that the complex function r(t) could follow. 



B. Zeroth Order action 



The zeroth order action is just 

S(a f , (/>*)= [ n> n^a^'dr = - r \2a^\a^'fdr (55) 
Jo Jo 

24/3 (V a) \i .24/3 

There are many different endpoints 77 which give the same a/, and many different con- 
tours to each of these end points. We note that these contours and end points are "physi- 
cally" different in that they correspond to different paths through the complex a space. Any 
one of these solutions can be graphed as a path through the complex a plane. That path 
is an invariant, and no coordinate transformation can alter that path. They thus represent 
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physically distinct paths. However, except for the parameters (3 , 7 and sign(af) differen- 
tiating the different solutions, the actual value of the action as a function of 0/ and aj is 
independent of which of the paths or of the endpoints are chosen. 

In the semi-classical wave function, the real part of the action determines the phase of 
the wave-function, while the imaginary part determines the amplitude. The two values of 7 
correspond to taking the complex conjugate of the iS while the two values of (3 correspond to 
two different values for the amplitude as well. These two possible values for the amplitude, 
given by the two possible values of /3, have been the subject of immense controversy. We will 
not take part in this controversy, but will keep (3 in the remaining equations. We do however 
note that the issue is not one of determining whether or not the "current" emitted by the 
a = "singularity" is purely outgoing or not. For both signs of (3 and 7 the system obeys 
the HH condition on the regularity of the wave-function at a = 0. Similarly, by taking an 
appropriate combination of the wave-functions with 7 = ±1 one can have a wave function 
with purely outgoing flux, purely in-going or a combination of the two. 

The usual path in r space, originally chosen in HH, is the path through r space beginning 
as a purely imaginary path from to —m/2v and then going parallel to the real axis to the 
endpoint [cosh" 1 {y a f) j v ^ — in/2u). a^(r) is real along the whole path. Along any other 
path, a has an imaginary part. We note that this path also corresponds to only one of the 
possible end values r . Other end points correspond to trajectories in the complex a plane 
which wrapping around a = more. In figure 2 we plot the trajectories in the complex a 
plane which correspond to the r trajectories of figure 1. 




Figure 2 

The paths in complex scale factor (a) space of solutions for the various complex time (r) 
paths of figure 1. 

The function, So(a,f,<f>f) is a Hamilton- Jacobi function for the e = system, in that it 
obeys the equation 



1 f dSo v 6o+ fj^y 2 +a3 , /o=0 . 



24a \ da J \2a 2 d(p / 

which is the Hamilton- Jacobi form of the constraint equation for e = 0. 



(57) 
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The first order action is just given by 



lS{a l ,h.Vo)= dS ° ia ^" V ° ) (58) 



This again solves the Hamilton- Jacobi function to first order 

+ a 3 = (59) 



1 dS (a<i))d5S(a,<i)) , 



12a da da 



C. First Order Solutions 

The first order equations for a and are required to find the second order action. Ac- 
tually, only the first order solution in is needed, but we will also solve the first order 
equations for a for future reference. The equations are 

- 12a<°> '5a' + 2V a^5a + a^ 2 (j> f = (60) 
( a (°) V)' = -a { ° )3 (61) 

with the boundary conditions that all of the variations are zero at r = 0. Because (f)^ is a 
constant, we also have that 50(t/) = 0. 
The solution for a is 

5a = (3~y- (it cos(iur) sm(iur)\ (62) 

while for 5(f) we get first that 



- [\ cos 3 (iz/r) - cos(iuT) + f ) 

= ^ ^ ' ^ ( 64 ) 

sin J (ivr) 

This expression has poles of order 3 at the points r = ^(2r + l)ir 

Note that at the end point, r = tq, the expression for 5a not only depends on the 
value of a/, but also on which particular end point r is chosen. This implies also that 
5t = ^ = — 5a(r ) /a^ (t )' will depend on which end point is chosen. Although neither 
of these enter into the expression for the second order action, we do expect them to have 
an effect on the third order action, giving an end-point dependence to the action to higher 
orders. 

Because 5(f) (to) = 0, we have 

/TO 
6<f>'df (65) 

v =-i/„ >,3<<f (66) 
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Because of the 1/a 3 term in the integrand for 5(f), its value will depend, among other things, 
on the way in which the path in r space wraps around the poles where a = 0. We get 



5<t>{r) 



2 

V 



1 + cos(iz/r) 



+ In 



1 + cos(iz/r) \ 1 

1 + cos(iuT ) J 1 + cos(iz/T ) 



(67) 



where the branches of the In function around the singularities where 1 + cos(iz/r) = are 
chosen so as to make the expression S4>(r(t)) continuous along the path r(t). 



D. Second Order Action 



Given the first order solutions, we can calculate the action to second order in e. 

d 2 



5 2 S = <f *S (a t , <P f , V ) - P a^\84>'fdr 

dVn JO 



We have 



Jo Jo ran (i 

i(3 



9z/ 6 



L 



IVT) 
cos(ii>To) ( (\ 



2\ 



- cos {ivr) — cos(iz/r) H — ) dr 



3/ 



z)(z 



1 + z 



dz 



(68) 

(69) 
(70) 



with z = cosiivr) which has poles of order 2 in the integrand at at z = — 1 or r = 
i(2m + 1)tt jv. Each time the integrand circles one of these poles, we accumulate a residue, 
and it is clear that the residue is the same at each of the poles in r space, since it is the 
same pole z = — 1 in z space. Thus we have 



Jo 



r ° a (0)3fXJj\2 • 



■z)(z + 2)' 
l + z 



96 

dz + 2nf3k—^ 



(71) 



where T is some fiducial path from the point r = to r = t , and k is the total net number 
of times that the actual path wraps around the various poles. (For an encircling of each 
of the poles in the counter clockwise direction, one gets the same contribution no matter 
which of the poles is chosen.) Note that the contribution from each of the poles to the 
action is real. They will therefore contribute only to the phase, and not the amplitude, of 
the wave-function. 

Assuming that va$ > 1, we can take the fiducial path as a straight line connecting r = 
to r = Tf. We finally get 



S 2 S = (3 



Si 



where Z = cos(iuT Q ) = ji\ 



2 d 2 S (a f , ( j )f ,Vol _ 7rjfc 192 



dV 2 
12 
Z+l 



12 In 



Z+l 



+ 17-12Z + Z C 



(72) 
(73) 
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We note that once again, the action does not depend on the endpoint of the integral-ie 
on which of the end values of Tq the integrand finishes on. 

We find that 5 2 S also obeys the second order Hamilton Jacobi function for the Einstein 
action. However it is important to note that to second order, S is a complex function 
of af,(f)f, with the imaginary part an actual function of these variables, instead of being 
a constant (as it is to zeroth order). Since the Hamilton Jacobi equation is a non-linear 
equation to second order in e, the real part of S will not be a Hamilton Jacobi function for 
the Einstein action. This will be important in the next section where we look at the classical 
paths which the wave-functions for the universe imply. 

III. CLASSICAL EVOLUTION 

The impression is often conveyed that the semi-classical path integral implies that the 
history that the universe follows is just the classical solution which is used to dominate the 
path integral. It is clear that this cannot be the correct answer. We have seen that the wave- 
function of the universe has contributions not just from real solutions to Einstein's equations, 
but also from complex solutions. These complex solutions clearly cannot represent a history 
of the universe, since physical things like the radius of the universe, the value of the scalar 
field, and their velocities in (proper) time are real quantities, not complex quantities. The 
semi-classical approximation to the path integral is a mathematical trick for evaluating the 
wave function, and does not in itself represent anything "real", and in particular does not 
represent a real history. However, once one has a wave function, one can ask the quantum 
mechanical question as to what types of history of the universe that wave function will 
represent. This is a very deep problem which does not have any satisfactory solution as yet. 
However, this paper is not the place for a detailed examination of the various possibilities. 
Instead we will take a naive approach. 

In ordinary quantum mechanics, the momentum operator conjugate to the position is 
just the imaginary derivative of the wave function. Although this operator really only makes 
sense if integrated over the configuration space, we will assume that some sense can be made 
of the momentum density, i/j(q)*i-^-ip(x). However this quantity is not real, while a physical 
momentum should be real. Furthermore, the amplitude squared of the wave function, |-^(g)| 2 
is the probability of finding the particle at the position q. The derivative of the amplitude 
of the wave function is thus not really something associated with the momentum of the 
particle. We can therefore take as the definition of the momentum at the point x, the 
quantity ip(q)i-j^ip(q), where ip(q) = iP(q)/\' 1 1'(q)\- Thus the momentum is just the derivative 
of the phase of the wave function. If the wave function is of the form e tS , then the momentum 
will be just the derivative of the real part of the action S. This argument is clearly highly 
suspect, and is at best only a crude approximation. However, it is the approximation which 
we shall use. Thus, the Real part of the action S will be used as a Hamilton Jacobi function 
for the evolution of the system as an approximation to the classical evolution implied by 
the quantum wave function. Again, we emphasize most strongly that this is at best a very 
crude attempt at deriving a classical evolution from the quantum system. In particular, 
in cases of significant quantum interference, this approximation will break down entirely, 
(for example, for a purely real wave-function, this gives the result that the momentum is 
everywhere zero). 
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The lowest order action is complex, but the imaginary terms do not depend on either a/ 
or <j)f. In the first order action, the imaginary part does depend on cf>f, but since the zeroth 
order action is entirely independent of <j)f, and the Hamilton Jacobi equation has ^ to the 
second order, this derivative of the imaginary part of the action will not contribute to the 
first order HJ equation. Thus, the real part of the first order action will again obey the HJ 
equation to first order. However, to second order, this is no longer true. The imaginary part 
of the first order action depends on <pf and thus comes in in quadratic order in the second 
order HJ equation. Ie, the real part of the second order action will not be a solution of the HJ 
equation, even though the full complex second order action is. Thus, the histories obtained 
by using only the real part of the action will not be solutions of the Einstein equations. 
Einstein's equations will have "quantum corrections" arising from the HH specification for 
the wave-function of the universe. 

In particular, if we use the real part of our second order solution, it obeys the HJ equation 
to second order in epsilon of 



+ 2 dS dS 2 S R 




(75) 

But in order for Sr to be a solution of the HJ equation for Einstein's equations to second 
order in e, the right hand side would have to be zero. Instead it is yz^s- We could subtract 
a term from the Einstein Hamiltonian constraint to find a new Hamiltonian which the real 
part would satisfy. We note that this term represents a potential which diverges as a goes 
to zero, a potential which will help prevent the universe from collapsing. 

H ^ = -A^ + A + a s VW - 2 SS-?- 3 (76) 

Note that the "quantum" correction falls off rapidly with a, which implies that once the 
universe is large, the impact of the imaginary parts of the action on the evolution of the 
universe become negligible. 

Let us examine the classical solutions obtained by using Sr as a HJ function for the 
solution. We are interested in the motion for v\df\ > 1, Furthermore, we expect Vq to enter 
into the complete action only in the combination V = V(<f)f) — Vq + e<pf. Thus, we define 




(77) 



and have 



24 

Sr « +PlyV 3 




+ 6i In 1 + llV + 12 7 y + jy 3 - 24tt k ) (78) 
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(Note that the term, 6i In could also be written as —12^arctan(y) which makes it clear 
that is is real). The equations relating the derivatives with respect to time of the dynamic 
variables and the HJ function are 

„ a = -12od = ^ (79) 
da 

where we have chosen the time to be the proper time in the space-time (ie, N=l). 
For large a, we have 



6 



1 + 6^ 



(81) 



and the evolution of ~ 2ey^rt, while a grows exponentially. These are just the "slow 
roll" results for the dynamics in an inflationary system, e makes only a small change in the 
results. 

What is more interesting is the behaviour of a and near the minimum, when V((f>)a 2 /6 
is of order unity. We first note that we can rewrite the expression for Sr as 

Sr - Ply ((1 + ^(y 2 ~ e 2 ^#) + e 2 0(y 5 ) + 0(e 3 ) (82) 

Ie, to second order in e and to fourth order in y this expression is identical to the previous 
expression for Sr, as can be readily verified by expanding this expression to second order in 
e. Since 

v 2 ~ = y fl2 " 1 " e2 ^n + 4 0q2 ~ Y (a " ar)(2 ° T) + 4^ + 0{{a ~ aT)2) (83) 



where 

a T = J^-(l + e 2 -^) (8J) 



6 M 2 2 



we can write the action as 



Defining £ = \a\ — ax, we have 



. 3 v 
2 \ 2 



(85) 



x. = -1204 « -12arf» ,8 7 |(1 + ^^(^ + 4# (*» 

^12^(1 + e^)(^ + 4# < 87 > 

^ = 4^^ 7 |M ( ^ + e | 0) i (88) 

«/»^(^ + ^)» (89) 
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where we have only kept terms to lowest order in £ and 0. Solving these equations we find 
that 

(^-{ + «f#)' = -3^(1 + ^)t (90) 

Since (-^y 2 ^ + e|0)^ was taken as being positive, we have that sign(t) = — Pj. Ie, the choice 
of (3 and 7 determines whether the solution starts or ends and that minimal value of a = ax- 
We also note that e only causes small (or order e 2 ) changes in the solution. 
We finally have 

1 

\a\ -or ~ {-a T + 0(e 2 ))t 2 (91) 
<t> » -e^ 2 (92) 

Note that the main effect of the HH boundary conditions has been to demand that is 
a maximum in time at the time when a is a minimum. One might have expected to be 
non-zero at this minimum, but it is not. This means that for these states, the value of 
cannot be used a clock near the turn around point, since the clock changes sign at point, 
although in principle could be. 

Of course, this turning point is also where one would expect to find maximal quan- 
tum interference between the various possible semi-classical paths, and so this non-analytic 
behaviour of the classical solution is suspect. However we note that the semi-classical wave- 
function, if we take an equal linear combination of terms with 7 = ±1, can be constantly 
interpreted not as a wave function which emerges out of the a = "singularity" (that would 
be consistent with taking only the solution with 7 = —/3), but rather with the boundary 
condition imposing the condition that the classical universe suffers a bounce. Ie, the HH 
conditions could be interpreted as enforcing a bouncing universe, rather than as describing 
the creation of the universe out of nothing. 

We note that we have neglected the contribution of the fluctuations to the path integral 
in all of the above analysis. One would expect that these too would be analytic in the 
path through r space which one chooses. However, these fluctuation terms could depend 
on which of the "end points" the path in r space ended at. In fact, because the first order 
solution 5a{r) is not periodic in imaginary r, one would expect the higher order actions to 
also depend on which of these end-points was chosen in defining the semi-classical action. 
Both the fluctuations and the end-point dependence could could substantially change the 
behavior of the wave function, especially near the "turning point" . 

IV. CONCLUSION 

We have solved the equations for the evolution of the universe in an inflationary approx- 
imation, where the slope of the potential driving the scalar field can be taken to be small. 
There exist an infinite number of complex solutions which can be said to dominate the path 
integral. Although some of these paths are analytically related to each other, others are not. 
The analytically related paths, which all have the same action, could, as others have done, 
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be regarded as equivalent in some sense. Since the manifold is real they are not related by 
a diffeomorphism, and the question as to whether or not they should be summed over or 
only one representative of the class should be used in defining the wave-function is to us still 
an open question. However, there exist singular points for the solutions in the complex r 
plane, which correspond to a going to zero. The path in the complex plane can wrap around 
these points, and the action depends on how the paths wrap around these singularities. To 
the order of approximation we examine, these singularities contribute to the real part of the 
path integral to second order in the slope of the potential. There exist an infinite number 
of these inequivalent analytic paths. Furthermore, there are also an infinite variety of ana- 
lytically inequivalent paths which depend on the choice of the end point of the integration 
in the complex r plane. The action, to second order, depends only on which of two classes 
the ends points fall into (as designated by the the parameter 7 above), although there are 
indications that to higher order the action will also depend on exactly which of these end 
points is chosen. 

If one assumes that the classical path that a universe would follow is best approximated by 
using the real part of the semi-classical complex action as the Hamilton Jacobi function, then 
one interpretation of the Hartle-Hawking wave function is not that the universe is created 
out of nothing, but rather that the HH condition forces the condition that the universe 
bounce rather than enter or leave the potential singularity. Of course it is precisely at the 
bounce point that the role of the quantum corrections (here taken to be the imaginary part 
of the semi-classical action) have their largest value, and the interpretation of the universe 
following a classical path is most problematic. 

We also note that our analysis demonstrates that , in the absence of non-homogeneous 
modes, the HH condition does not select the sign of the real part of the action (the parameter 
f3 above) and thus this analysis does not have anything to say to the controversy over this 
sign. This model however seems simple enough that it could also be used to address the 
issue of non-homogeneous fluctuations, at least to some low order of approximation, which 
could cast some light on the controversy. However, we believe that this would require a 
deeper understanding than we at least have of the role that such complex solutions really 
play in the evaluation of a path integral. 

What has not been done here is to calculate the contribution of the " determinant" , (the 
integration over non-classical fluctuations), to the wave function. It is possible that this 
would also lift the degeneracy over the various possible paths. This paper also does not 
address the broader question of the how one should use such complex solutions in evaluating 
the path integral. However, even in the context of ordinary quantum mechanics with a single 
degree of freedom, this question of how complex classical paths are to be used in general in 
evaluating the propagator is still one requiring more understanding. 
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